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ABSTRACT 

The analysis of the CoRoT space mission data was performed aiming to test a method that 
selects, among the several light curves observed, the transiting systems that likely host a low- 
mass star orbiting the main target. The method identifies stellar companions by fitting a model 
to the observed transits. Applying this model, that uses equations like Kepler's third law and 
an empirical mass-radius relation, it is possible to estimate the mass and radius of the primary 
and secondary objects as well as the semimajor axis and inclination angle of the orbit. We 
focus on how the method can be used in the characterisation of transiting systems having a 
low-mass stellar companion with no need to be monitored with radial-velocity measurements 
or ground-based photometric observations. The model, which provides a good estimate of the 
system parameters, is also useful as a complementary approach to select possible planetary 
candidates. A list of confirmed binaries together with our estimate of their parameters are 
presented. The characterisation of the first twelve detected CoRoT exoplanetary systems was 
also performed and agrees very well with the results of their respective announcement papers. 
The comparison with confirmed systems validates our method, specially when the radius of the 
secondary companion is smaller than 1.5 Rjup, in the case of planets, or larger than 2 Rjup, in 
the case of low-mass stars. Intermediate situations are not conclusive. 

Subject headings: planetary systems — techniques: photometry — techniques: transit mod- 
elling 



' The CoRoT space mission, launched on 2006 December 27, was developed and is operated by the CNES, with participation 
of the Science Programs of ESA, ESA's RSSD, Austria, Belgium, Brazil, Germany and Spain. 
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1. Introduction 



The analysis of transiting systems based on light curve modelling combined with ground-based follow- 
up by means of radial-velocity measurements have shown their exceptional importance in the characteri- 
sation of extrasolar systems. Together, they provide the determination of several physical and orbital pa- 
rameters that cannot normally be obtained for non- transiting systems, such as the mass and radius of the 
secondary companion, and thus its density. The possibility of a detailed study of extrasolar transiting light 
curves in a large number of targets has been co nsiderably improved by recent photometric space missions, 
such as CoRoT dBaghn et alJboOa) and Kepler dBorucki et al.ll2010l) . 



At the time of this writing the CoRoT space telescope has collected, since its launch, the light curves 
of more than a hundred thousand stars through 13 observational runs. This huge number of stars is a strong 
motivation to develop tools to efficiently treat the released data, specially considering that first of all the data 
need to be cleaned from long-term variations, short-term oscillations, outliers, discontinuities, and others. 
Regarding the characterisation of extrasolar systems, an extra challenge is to distinguish transits caused by 
planetary companions from those related to the presence of a low-mass star in a binary system, particularly 
given the time-consuming ground-based observations that are normally used as a complementary approach. 



Carpano et al.l (120091) published their results concerning the analysis of CoRoT light curves observed 
during the initial run of the mission, named IRaOlE in whi ch they presented a l ist of 50 planetary candidates 
together with a hst o f 145 eclipsing binary candidates. Moutou et al. fcoo ob complemented the work of 
Carpano et al.l (|2009|) wit h additional follow -up observations, also in the context of the initial run. In other 
two papers, published by IPeeg et all (|2009h and lCabrera et al.l (|2009h . the authors presented a list of targets 
for which ground-based follow-up was conducted, helping as a complementary approach in the classification 
of the candidates. The lists released by these papers are the result of a huge effort of several working groups, 
and shows how difficult it is to characterise systems, planetary candidates or not, among thousands of targets. 



In the work of ISilva & Crua (l2006h . the authors proposed a method based on the fit of light curves with 
transits that can be used in the characterisation of transiting systems. This method provides the determination 
of some parameters of the system, such as the orbital inclination angle, the semimajor axis, and the mass 
and radius of the primary and secondary objects. In the present work we have used an updated version of the 
same method, applied to a list of transiting light curves from publicly available runs of the CoRoT mission. 
The purpose of our analysis is to show that the method is useful to identify transits most likely caused by a 
binary configuration, without making use of any time-consuming effort to conduct ground-based follow-up, 
like radial velocity measurements or photometric observations. 



ISilva & Cruzl(|2006f) also tested their model applied to the transiting exoplanetary systems known at that 
time, and their results are in good agreement with the published parameters. Therefore, here we analysed 
the first confirmed CoRoT planetary systems as well, from CoRoT-1 through 12, comparing our estimates 
to those in their respective announcement papers. 



^ IR: Initial Run; LR: Long Run; SR: Short Run; c and a represent the direction of tiie galactic and tiie anti-galactic centre, 
respectively. 
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Fig. 1. — Example of polynomial function fit perfomied to correct the light curves from eventual disconti- 
nuities. Only points within twice cr and located outside transits (dots) are taken into account when searching 
for the best fit (solid lines). Transit positions are represented by crosses. This figure represents CoRoT-1, 
which folded light curve is shown in Fig. |6] after the discontinuity corrections were made. 



Section |2] presents the data reduction, showing the corrections needed to apply to the light curves (in 
the format delivered to the scientific community) before modelling the transit shape. Next, Sect.|3]describes 
the method and its usefulness in the classification of binary system candidates. In Sect. HI the results are 
presented and discussed, which includes our parameter estimate for a list of binary systems first identified 
by other works as well as the characterisation of confirmed CoRoT exoplanetary systems. Finally, Sect. |5] 
presents the final remarks and conclusions. 



2. Sample data and reduction 



The light curves analysed using our method are part of the three-colour band data from publicly avail- 
able runs observed by the CoRoT mission (a few cases of monochromatic light curves were also included). 
Before doing any kind of fit to model the observed transits, the light curves have to be cleaned from any 
systematic noise th at may still remain in the data delivered to the scientific community, which format is 
called the N2 level (jBaudin et al.ll2006r) . The systematic noise normally seen is related to: /) discontinuities 
produced by hot pixels; //) outliers, whose sources are diverse; and ///) short-term oscillations related to the 
Co RoT orbital frequency (103 rn in) a nd its harmonics. For details on the CoRoT satellite and its orbit see 
e.g. iBoisnard & Auvergnel (120061) and lAuvergne et al.l(l2009h . In addition to the systematic noise, short-term 
oscillations, intrinsic to some type of stars, were also observed in some light curves and removed. 



We developed a code to correct the light curves from discontinuities and outliers. The code computes 
the first derivate of the data points in order to locate the discontinuities, then it fits polynomial functions to 
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Fig. 2. — Phase-folded light curve with short-term oscillations (in this case due to the CoRoT orbit) that we 
have to correct before modelling the transit shape. The same system is shown in Fig.[5]after these oscillations 
have been removed. 

points between them, and finally the light curve is normalised according to such functions. When searching 
for the best fit, the code does not consider: /) data points located outside twice cr, which value is estimated in 
a specific portion of the light curve; and //) data points located where transits happen. We note that a visual 
inspection is always conducted to avoid any abrupt behaviour of the fitted function, specially close to transit 
times. Figured] shows an example of this polynomial function fit. 

Short-term oscillations may also affect the search for the model that best fits the observed transit. 
Figure |2] shows an example of a transiting light curve when this kind of oscillation is present. The same 
object is shown in Fig. |5] after being corrected from the most prominent harmonics. Oscillations due to the 
intrinsic variability of some stars were also properly removed when needed. 



3. Light curve fit 



The method developed by Silva & Cruz ( 200 6^ was used here to search for the best model parameters 
that fit the observed light curve transits. The model considers an opaque disc that simulates the secondary 
object passing across the stellar disc. We assumed a quadratic function to describe the limb-darkening of 
the disc of the prim ary object, which is based on the star HD 209458 (ui = 0.2925 and U2 = 0.3475, from 



Brown et al.ll200lh . An exception is the system 0100773735 (LRcOl, Fig. [5]l, for which ui and U2 were 



assumed to be half those of the star HD 209458. In the case in which the secondary companion is a low- 
mass star, it will not be an opaque disc. However, this will not considerably change the results, since its 
flux contribution is small compared to the main star (a 0.3 M© star orbiting a solar-type star contributes less 
than 2% to the total flux). Cases in which the secondary is as bright as the primary were not selected in our 
analysis given the large transit depth that would be observed in the light curve. 



The orbital period (P) of the companion is a known parameter, obtained directly from the light curve, 
whereas the following three variables are the result of the best fit: the radii ratio between secondary and 
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Fig. 3. — Example ofx^ minimisation for the system 0100773735 (LRcOl, Fig.[5]l in the search for values 
of Rp, dp, and / of the model that best fit the observed transit. 

primary objects (Rp = R2/R\), the semimajor axis of the secondary orbit in units of the primary radius 
(flp = a jR]), and the orbit al inclination angle (/). The search for the best fit is conducted with the AMOEBA 
routine ( Press et al which performs a multidimensional chi-square ix^) minimisation of the function 

f{Rp, Up, i) describing the transit profile. 

The first set of parameters normally used as input are: Rp - 0.1, / - 85°, and Qp = a/Ri calculated 
according to the Kepler's third law (Equation |3]l for M\ + M2 ~ 1 M© and R\ ~ I R©. In order to avoid 
premature convergences on local minimums and explore the entire parameter space, the best solution is 
found after running the routine several times, replacing the input parameters by new values within a chosen 
range (e.g. ARp = ±0.05, Aop = ±3.0 and A/ = ±5°). Throughout the execution of the process, a visual 
inspection of the fit is carried out after each possible solution is achieved. Figure |3] shows an example of 
minimisation for the system 0100773735 (LRcOl). 



3.1. Estimate of mass and radius 



The method used to estimate the mass and radius of the primary and secondary companion was up- 
dated using new mass-radius relations based on more recent discoveries, specially the one fitted to known 
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Fig. 4. — Empirical mass-radius relations used by our method. The long-dashed line is a function, given by 
Equation [H fitted to stars (open circles) whereas the solid line represents a function, given by Equation |2] 
fitted to known exoplanets (open squares) and the planets of the Solar System (filled squares). Crosses 
represent exoplanets and low-mass objects that were not used in the fit (see text). The values of mass and 
radius were taken from the Extrasolar Planets Encyclopaedia^. 

exoplanets (systems presently listed in the Extrasolar Planets Encyclopaedic^] were used). These relations 
are: 

R - AM^ (1) 

R = CM + D (2) 

where A = 0.013938, B = 3.867270, C = 0.410162, and D = 0.065934. The function given by Equation^ 
was fitted to stars whereas Equation |2]represents a function fitted to known exoplanets and the planets of the 
solar system. Both functions are plotted in Fig. |4l Exoplanets having radius > 1 Rjup or mass > 2 Mjup and 
low-mass objects having radius < 2 Rjup were not used in the search for the best fit. These objects belong to 
an ambiguous region of the mass-radius diagram, in which a change in mass does not necessarily imply in a 
change in radius. 

After applying our method to the observed transits, we are able to estimate the mass and radius of both 
the primary (Mi, /?i) and secondary (M2, R2) objects and the orbital semimajor axis a (given in astronomical 
units) using the following relations: 

= —r{Mi + M2) (3) 



47^2 

R2 



^ |http://exoplanet.eu| 
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a 



(5) 




R2 ^ AM2 



B 



(6) 



Ri =AMi", 



R2 = CM2 + D 



(V) 



where Equation |3]is the Kepler's third law, Equations |4] and |5] are obtained by the transit best fit, and Equa- 
tions |6] and |7] are the empirical mass-radius relations used depending on whether the secondary companion is 
a star or a planet. Therefore, two sets of five relations aie numerically solved and two sets of the parameters 
Ml, Ri, M2, R2, and a are computed. If R2 > 2 Rjup then we consider the system to be a binary candidate 
and the parameters yielded by Equations |3] S] |5] and |6]are used. On the other hand, if /?2 < 2 Rjup, then the 
parameters yielded by Equations |3] HI |5] and Hare used. 

The numerical calculation of M\, R\, M2, R2, and a proceeds as follows: first, one of the parameters 
is fixed (e.g. Mi = I Mq) and the others are calculated; then, the ratio a/Ri is compared to the value 
of Up provided by the transit best fit; the fixed parameter is iteratively incremented (or decremented) until 
the difl'erence between a/Ri and ap is as smaller as one wishes. When this condition is satisfied, the five 
parameters are finally estimated. 

Concerning the mass of the secondary objects (M2), there are three possibilities: /) if the companion 
has R2 > 2 Rjup, then it will probably be a low-mass star and its mass can be easily estimated by the mass- 
radius relation fitted to that region of the diagram; //) if /?2 < 1 Rjup> then we can estimate a value for M2 
using the other mass-radius relation; however, since the mass-radius relation for planets depends on their 
chemical composition, it is not possible to estimate an accurate mass value for the secondary in this radius 
regime if its composition is unknown; ///) finally, if 1 < R2 < 2 Rjup, then the mass can not be univocally 
determined since, for one given radius, objects in this radius regime may have masses ranging from about 
1 Mjup to about 100 Mjup (brown dwarfs). Based on these arguments, we only include in our discussion a 
mass estimate for secondary companions with R2>2 Rjup. Nevertheless, we note that our method provides 
a good estimate for the radius of not only primary and secondary stars, but also of planetary companions 
(even for R2 < 2 Rjup). Indeed, this is the case of the confirmed CoRoT systems, which results are presented 
and discussed in Sect. ID 



The results presented here are for the parameters a, i, Mi,Ri, M2, and R2, where a is given in astronom- 
ical units. Ml and M2 in solar masses, Ri in solar radii, and R2 in Jupiter radii. To estimate the uncertainty 
in these parameters, the procedure is as follows: 



3.2. Uncertainties in the system parameters 



1 ) first, the standard deviation (cr) is computed in a region outside a transit for the model that best fits 
the light curve; in a region inside a transit, a possible pre sence of spots on the surface of the primary 
star would cause variations in the light curve (|Silvall2003|) and cr would be miscalculated; 
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2) by changing, at a given step, the three basic parameters yielded from the best fit to the transit {Rp, Up, 
and /), we obtain a new light curve; 

3) next, these parameters are changed until the new light curve starts deviating significantly from the 
best one; that is when the difference between the light curves is more than 1 cr for at least 10% of the 
transit; 

4) finally, these new basic parameters are used to re-estimate M\, M2, Ri, and/?2 following the procedure 
described in Sect. 13.11 the difference between both new and best estimates provides the uncertainties 
in the parameters. 



Our mass-radius relations (Equations [T]and|2]i do not depend directly on the orbital inclination angle. 
However, since a change in this parameter leads to a change in both depth and duration of the transit, its 
uncertainty can be included in the uncertainty determination of the other two basic parameters. To do so, 
first the uncertainty in the inclination angle is estimated as described above (in steps 2 and 3 only / is 
changed). Then, the uncertainties in the other two parameters (each one in turn) are estimated considering 
the inclination angle changed by its uncertainty. At the end, the uncertainties in M\, Ri, M2, and R2 take into 
account those in Rp, ap, and /. The uncertainty in a includes those in both ap and Ri given the conversion of 
the semimajor axis from units of stellar radius to astronomical units. 

To reduce noise and clarify visualisation, a smoothing function was applied to the light curves before 
searching for the model that best fits the transit. In addition, for some systems only the best quality channels 
were used in the analysis. These procedures result in a more accurate estimate of the system parameters and, 
consequently, in smaller uncertainties. 



4. Results and discussion 

4.1. Characterisation of binary system candidates 

Tableland Fig. [5] show a few examples of our method applied to light curves of systems classified as 
eclipsing binaries or targets having faint background eclipsing binary stars. 



According to iSilva & Cruzl (120061) . the present method should consider as non-planetary candidates 
only the systems for which the radius of the secondary companion is larger than 1.5 Jupiter radii. After that 
publication, some exoplanets with radius between 1.5 and 2 Rjup were discovered. Thus, in this work, we 
chose a more conservative value of /?2 - 2 Rjup for the lower limit of binary candidates. 

The systems that our method classify as binary candidates actually represent classes of systems with 
different configurations of the components. There are several sources of false alarms that can mimic the 
transit of a planet in front of the main target. The most common are eclipsing binary systems with grazing 
transits, targets having a background eclipsing binary system, or even triple systems. The present method 
does not identify the exact configuration of the system, but provides the information that the observed transits 
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Table 1 : Physical and orbital parameters for CoRoT systems analysed in this work. The uncertainties in a, i, 
Mi,Ri, M2, and R2 were estimated as described in Sect. 13.21 These targets represent confirmed or probable 
binary systems and they are all commented in Sect. 14.11 
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a 


/ 


M, 




M2 


i?2 




CoRoT ID 


Run 




[AU] 








Ref. 


[days] 


[deg] 


[Mo] 


[Re] 


[Mg] 


[Rjup] 


0102787048 


IRaOl 


7.896 


0.094 ± 0.008 


85.1 ±0.2 


1.63 ±0.08 


2.1 ±0.2 


0.16 ±0.07 


2.00 ±0.31 


[2, 3, 4] 


0102811578 


IRaOl 


1.66882 


0.0339 ± 0.0008 


77.1 ±0.2 


1.57 ±0.02 


1.95 ±0.04 


0.32 ± 0.02 


3.09 ±0.17 


[2] 


0102815260 


IRaOl 


3.587 


0.057 ± 0.006 


> 87.0 


1.71 ±0.09 


2.2 ±0.2 


0.19 ±0.08 


2.14 ±0.31 


[2, 4] 


0102855534 


IRaOl 


21.72 


0.218 ±0.014 


86.9 ± 0.2 


2.48 ± 0.07 


4.2 ± 0.2 


0.50 ± 0.05 


4.49 ± 0.42 


[2, 4] 


0100773735 


LRcOl 


4.974 


0.074 ± 0.004 


82.5 ± 0.2 


1.98 ±0.06 


2.8 ±0.1 


0.23 ± 0.04 


2.42 ± 0.24 


[1] 


0100885002 


LRcOl 


11.8054 


0.136 ±0.007 


85.1 ±0.1 


2.02 ± 0.02 


2.9 ±0.1 


0.42 ± 0.04 


3.80 ±0.31 


[1] 


0101482707 


LRcOl 


39.89 


0.270 ± 0.004 


88.5 ±0.1 


1.45 ±0.01 


1.73 ±0.02 


0.19 ±0.01 


2.18 ±0.09 


[1] 


0101095286 


LRcOl 


5.053 


0.088 ± 0.005 


70.6 ± 0.4 


3.2 ±0.1 


6.7 ± 0.3 


0.43 ±0.12 


3.91 ±0.92 


[1,3] 


0101434308 


LRcOl 


79.95 


0.41 ± 0.02 


> 89.7 


1.30 ±0.05 


1.46 ±0.08 


0.18 ±0.06 


2.10 ±0.17 


[1] 


0211660858 


SRcOl 


8.825 


0.114 ±0.008 


86.3 ±0.3 


2.12 ±0.07 


3.2 ± 0.2 


0.46 ± 0.05 


4.10 ±0.38 




0211654447 


SRcOl 


4.751 


0.090 ± 0.005 


70.9 ± 0.4 


3.3 ±0.1 


7.4 ±0.3 


1.04 ±0.10 


9.90 ± 1.23 




0102755837 


LRaOl 


27.955 


0.28 ± 0.02 


83.8 ± 0.2 


3.4 ±0.1 


7.6 ± 0.4 


0.48 ± 0.07 


4.28 ± 0.61 





[l]|Cabreraet"Z|j2009t); [2]|Carpano e7ai](l2009h; [3]|Deeg et al]j2009t); [4]|Moutou etal] |2009|). 



are likely not caused by a planet. Listed below are our comments for each case. The window ID of the run 
(e.g. El -028 8) is also shown as a complementary identifier. 



IRaOl - 0102787048 (El-0288) 



This system was first class ified by ICarpano et al.l (120091) as a planetary transit candidate. However, af- 
ter follow-up observations, iMoutou et al.l (|2009b confirmed that the transit is originated by a background 
eclipsing binary and then diluted by the main target. They estimated a mass-ratio of 0. 15 between secondary 
and primary components of the binary system. Using the results of our method we have a value of 0.10 + 
0.04 for the same ratio, which is consistent with a diluted transit of a more massive object. The fact that no 
transit is observed in the red channel also helped to label a non-planetary nature for this target. 



IRaOl - 0102811578 (E2-0416) 



System classified by ICarpano et all (l2009h as an eclipsing binary. No information concerning follow-up 
observations has been published. Our estimate of mass and radius for the secondary companion confirms its 
binary nature. 



IRaOl - 0102815260 (E2-2430) 



Also i n the list of planetary transit candidates of lCarpano et all (120091) . but afterwards classified by lMoutou et al. 
(120091) as a binary system according to radial- velocity observations. Our estimate for the mass-ratio is 0. 11 
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± 0.04, slightly smaller than the value of 0.17 published by iMoutou et all (120091) . but still consistent with a 
stellar companion. 



IRaOl - 0102855534 (E2-1736) 



Follow-up observations conducted by iMoutou et al.l (120091) indicate th at an eclipsing binary is the main 
target, changing the planetary nature of the secondary first suggested by ICarpano et al.l (120091) . The radial- 
velocity measurements indicate a mass-ratio of 0.2, which agrees with our estimate of 0.19 ± 0.01. Indeed, 
the mass and radius that we show in Table [T] clearly indicate the stellar nature of the secondary. 



LRcOl - 0100773735 (E2-1245) 



Based on radial-velocity observations, ICabrera et al.l (120091) suggest that this is a spectroscopic binary or 
multiple system. Our mass and radius estimate confirms their conclusion of a non-planetary object causing 
the observed transits. 



LRcOl - 0100885002 (E2-4653) 



This target is listed in ICabrera et al.l (|2009[) as an eclipsing binary, which is in agreement with the mass 
and radius estimated using our method. Transits are observed only in the blue channel, contributing to the 
classification of this system as non-planetary. 



LRcOl - 0101482707 (El-2837) 



Also classified by ICabrera et al.l (|2009l) as a binary system after radial- velocity observations and confirmed 
by our method. 



LRcOl - 0101095286 (El-2376) 



The follow-up of this target with photometric observations unveiled its binary nature (jCabrera et al.ll2009|) . 
which can be clearly confirmed by our results of mass and radius for the secondary object. 



LRcOl - 0101434308 (El-3425) 



ICabrera et al.l(|2009 l). without doing any follow-up observations, concluded that this is a binary system given 
the fact that the transit is predominantly seen in the blue channel. 

In Tabled] are also listed our results for the targets 0211660858 (SRcOl, E2-0369), 0211654447 
(SRcOl, El-1165), and 0102755837 (LRaOl, E2-2249). The CoRoT team has not yet pubUshed the re- 
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suits of their analysis for these runs. Nevertheless, they are shown here because the secondary mass and 
radius estimated by our method clearly classify these targets as non-planetary systems. 

As one can see in Table [T] for some targets the estimated radius of the secondary companion is close 
to the limit of 2 Rjup used to distinguish binaries from possible planetary systems. They were classified 
as probable binary systems considering that, at present, all detected exoplanets for which the radius was 
derived are as large as 1.8 Rjup or smaller. 

Among the nine targets listed in Table [Uthat were analysed and published by the CoRoT team, six were 
classified as binary systems only after ground-based follow-up, either by radial-velocity measurements, or by 
photometric observations, or both. This is an indication that the method normally used by different CoRoT 
working groups is perhaps not good enough to provide a pre-classification of the targets. Our method would 
exclude the binary targets by itself, and no time-consuming follow-up observations would be required for 
most cases. 



4.2. Characterisation of exoplanetary systems 

Table |2] presents our parameter estimates for confirmed CoRoT exoplanetary systems, from CoRoT- 

1 through 12, together with the results presented in the announcement papers. Here the method uses the 
information provided by the mass-radius relation fitted to known transiting exoplanets and to the planets of 
the solar system (Equation |2] Fig. ID. Figure |6]shows the light curves of the first six CoRoT systems plotted 
with our best fit. 

The results of a detailed comparison between our estimates and those of the announcement papers are 
shown in Table|3] This table lists the coefficients of weighted linear regressions obtained for a, i, M\,R2, and 
R2, where the weights 1 /cr^ were used for cr representing the errors that we estimated for these parameters. 
At first, we computed the coefficients using the first twelve CoRoT systems and no systematic difference 
was found within 2cr. Within Icr, however, the agreement is not so good, specially for M\ and R\. This is 
most due to discrepancies in the comparison of CoRoT-8, as one can observe in Table |2] The data points in 
the light curve of this system has a dispersion of 0.017, which is much larger than those of the other CoRoT 
systems: -0.007 or smaller. Indeed, if the parameters of CoRoT-8 are not included in the regressions, 
the agreement between our results and those of the announcement papers is much better, standing within Icr 
mostly. An extra caution when dealing with noisy light curves is, therefore, recommended. Light curves with 
small-depth or long-period transits (such as CoRoT-7 and 9, respectively) also produce larger uncertainties 
in the parameters and should be carefully analysed. 

In Table 121 all planets detected in the CoRoT field have /?2 < 1.5 Rjup, also including our estimate for 
CoRoT- 1 if we take its uncertainty into account. This result supports our suggestion that the present method 
can be used to characterise not only systems for which the radius of the secondary companion is larger than 

2 Rjup, but also those having R2 smaller than 1.5 Rjup, helping as a complementary approach in the search 
of promising candidates for radial- velocity follow-up. The intermediate situations, when R2 stands between 
1.5 and 2 Rjup, are not conclusive. 
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Table 2: Parameters of the first twelve confirmed CoRoT planetary systems compared to the results of this 
work. The period values are from the respective announcement papers. The uncertainties in a, i, M\, R\, 
and /?2 were estimated as described in Sect. 







P 


a 




Ml 


R\ 


R2 




CoRoT ID 


Run 


[days] 


[AU] 


[deg] 




Ref. 


[Mo] 


[Ro] 


[Rjup] 



0102890318 

(CoRoT-1) 

0101206560 

(CoRoT-2) 

0101368192 

(CoRoT-3) 

0102912369 

(CoRoT-4) 

0102764809 

(CoRoT-5) 

0106017681 

(CoRoT-6) 

0102708694 

(CoRoT-7) 

0101086161 

(CoRoT-8) 

0105891283 

(CoRoT-9) 

0100725706 

(CoRoT-10) 

0105833549 

(CoRoT-11) 

0102671819 

(CoRoT-12) 



IRaOl 
LRcOl 
LRcOl 
IRaOl 
LRaOl 
LRc02 
LRaOl 
LRcOl 
LRc02 
LRcOl 
LRc02 
LRaOl 



1.5089557 (64) 
1.7429964(17) 
4.256800 (5) 
9.20205 (37) 
4.0378962 (19) 
8.886593 (4) 
0.853585 (24) 
6.21229 (3) 
95.2738 (14) 
13.2406 (2) 
2.994330(11) 
2.828042 (13) 



0.027 ± 0.002 


85.2 ±0.5 


1.11 ± 0.04 


1.18 ±0.06 


1.59 ±0.13 


[0] 


0.0254 ± 0.0004 


85.1 ± 0.5 


0.95 ±0.15 


1.11 ±0.05 


1.49 ±0.08 


[1] 


0.028 ± 0.002 


87.5 ± 0.4 


0.90 ± 0.03 


0.90 ± 0.04 


1.38 ±0.10 


[0] 


0.0281 ±0.0009 


87.8 ±0.1 


0.97 ± 0.06 


0.902 ±0.018 


1.465 ±0.029 


[2] 


0.056 ± 0.005 


86.3 ± 0.4 


1.33 ± 0.07 


1.52 ±0.11 


0.98 ± 0.13 


[0] 


0.057 ± 0.003 


85.9 ±0.8 


1.37 ±0.09 


1.56 ±0.09 


1.01 ± 0.07 


[3] 


0.087 ± 0.006 


> 89.3 


1.03 ± 0.04 


1.07 ± 0.06 


1.07 ± 0.10 


[0] 


0.090 ± 0.001 


> 89.915 


1.16 ±0.03 


1.17 ±0.03 


1.19 ±0.06 


[4] 


0.053 ± 0.003 


85.0 ± 0.2 


1.19 ±0.04 


1.29 ±0.06 


1.37 ± 0.14 


[0] 


0.04947 ± 0.00029 


86 ± 1 


1.00 ± 0.02 


1.19 ±0.04 


1.388 ± 0.047 


[5] 


0.083 ± 0.006 


89.4 ± 0.4 


0.96 ± 0.05 


0.98 ± 0.06 


1.06 ±0.10 


[0] 


0.0855 ± 0.0015 


89.1 ± 0.3 


1.05 ± 0.05 


1.025 ± 0.026 


1.166 ±0.035 


[6] 


0.018 ±0.005 


78.2 ± 1.5 


0.98 ± 0.17 


1.01 ± 0.23 


0.15 ± 0.07 


[0] 


0.01720 ± 0.00029 


80.1 ± 0.3 


0.93 ± 0.03 


0.87 ± 0.04 


0.150 ±0.008 


[7] 


0.067 ± 0.004 


86.7 ±0.1 


1.06 ±0.03 


1.11 ± 0.05 


0.88 ± 0.09 


[0] 


0.063 ± 0.001 


88.4 ±0.1 


0.88 ± 0.04 


0.77 ± 0.02 


0.57 ± 0.02 


[8] 


0.39 ± 0.02 


> 89.87 


0.85 ± 0.03 


0.84 ± 0.04 


0.93 ± 0.08 


[0] 


0.407 ± 0.005 


> 89.95 


0.99 ± 0.04 


0.94 ± 0.04 


1.05 ± 0.04 


[9] 


0.108 ±0.005 


88.0 ± 0.1 


0.94 ± 0.03 


0.95 ± 0.04 


1.15 ± 0.10 


[0] 


0.1055 ± 0.0021 


88.6 ±0.2 


0.89 ± 0.05 


0.79 ± 0.05 


0.97 ± 0.07 


[10] 


0.044 ± 0.002 


83.2 ±0.2 


1.24 ± 0.04 


1.37 ±0.06 


1.32 ±0.13 


[0] 


0.044 ± 0.005 


83.17 ± 0.15 


1.27 ±0.05 


1.37 ±0.03 


1.43 ±0.03 


[11] 


0.040 ± 0.003 


85.7 ±0.2 


1.02 ±0.05 


1.06 ±0.07 


1.33 ± 0.15 


[0] 


0.0402 ± 0.0009 


85.5 ±0.8 


1.08 ±0.08 


1.12 ±0.10 


1.44 ±0.13 


[12] 



rOl this work: [lllBarae et al.l fc008l): [21 


Alonso et a 


.I(l2008b: [31 


Deleuil et aljj2008h: [4llAi2rain et al.ll2O08l); 151 


RaueretalJll2009l) 


[6] Fridlund et al. (2010); [7] Leger et al. 


i2009|); [81 


Borde et alj 


2010l): [9llDee2 et alj J2010l): [lOllBonomo et alj 


20101); 



5. Conclusions 

We have presented a method that provides a good estimate of some physical and orbital parameters of 
a transiting system, such as the mass and radius of the secondary companion. Applied to transiting light 
curves, the method will exclude cases most probably related to low-mass stars in a binary system, instead 
of a planet. In other words, our method is able to exclude systems that at first may be considered as good 
planetary candidates but that afterwards would have their binary nature unveiled, without making use of 
time-consuming ground-based measurements normally conducted to complement the observations. 

We note that the method does not, by itself, determine the real nature of the secondary object (whether 
it is a binary companion or not). Instead, it identifies and characterises probable candidates for binary 
systems, which will help to reduce the huge number of targets initially available and to create a list of 
priority stars, still candidates for planetary systems, to be monitored with radial-velocity measurements. We 
do not discard other methods, however, which can be used to complement our approach. 
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Table 3: Coefficients of the linear regression obtained in the comparison of the parameters listed in Table |2] 
with CoRoT-8 included or not (see Sect. 4.2). The angular and linear coefficients (a, /?), the dispersion 
around the fit {nns), and the correlation coefficient (cc) are shown for both scenarios. 

a p nns cc 

For the first twelve CoRoT systems: 

a 0.98 + 0.02 0.0016 + 0.0010 0.0017 1.00 

0.91 ±0.12 7 ± 10 0.68 0.97 

Ml 0.55 ±0.27 0.47 ± 0.28 0.11 0.68 

Ri 0.66 ±0.19 0.38 ±0.19 0.12 0.84 

R2 0.88 ±0.09 0.11 ± 0.10 0.13 0.94 

CoRoT-8 not included: 

a 0.97 ±0.02 0.0017 ± 0.0009 0.0015 1.00 

0.97 ±0.06 2± 5 0.39 0.98 

Ml 0.72 ±0.28 0.26 ± 0.29 0.10 0.73 

R, 0.84 ±0.16 0.17 ± 0.16 0.09 0.91 

R2 0.95 ± 0.07 0.02 ± 0.07 0.09 0.97 



The method was also applied to twelve CoRoT targets confirmed as planetary systems, showing that the 
estimated radii of the secondary companions (as well as other orbital parameters) are in very good agreement 
with the results published by the respective announcement papers. This means that our model can also be 
used in the characterisation of possible exoplanetary systems, specially when R2 is smaller than or of the 
order of 1.5 Rjup. No conclusions could be drawn concerning the radius range 1.5 < /?2 < 2 Rjup- 

Our model is useful not only to be applied to CoRoT light curves that have been or will be released 
by the mission, but also to data of other present or future missions based on photometric observations of 
transiting systems that involve a large sample of targets (such as the Kepler mission). 

We thank the financial support fromFundagao de Amparo a Pesquisa do Estado de Sao Paulo (FAPESP) 
in the form of a grant (2006/50654-3) and a fellowship (2008/03855-9). We also thanks the Instituto Na- 
cional de Pesquisas Espaciais (INPE) for its support. 
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Fig. 5. — Examples of binary systems analysed in the present work. Here and in the following figure the 
light curves were smoothed to reduce noise and clarify visualisation. The index c, when added to the CoRoT 
ID, indicates that the light curve was corrected from short-period oscillations. 
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Fig. 6. — The first six CoRoT exoplanetary systems: CoRoT-1 tiirougli CoRoT-6, respectively, from left to 
right and from up to bottom. 
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